Surface spin waves in superconducting and insulating ferromagnets 
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Surface magnetization waves are studied on a semi-infinite magnetic medium in the perpendicular 
geometry. Both superconducting and insulating ferromagnets are considered. Exchange and dipole 
energies are taken into account, as well as retardation effects. At large wave vectors, the spectrum 
for a superconductor and insulator is the same, though for the former the branch is terminated much 
earlier than for the latter due to excitation of plasmons. At small wave vectors, the surface wave is 
more robust in the superconductor since it is separated from the bulk continuum by a finite gap. 
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Studies of spin waves began long ago and proved to be 
a powerful method of investigation of magnetically or- 
dered materials Q . These studies can be divided into two 
categories. In the first category, problems are included 
regarding spin waves in the bulk of an infinite medium, 
where long-range geometric effects do not arise, while the 
spin-exchange stiffness is taken into account. The second 
category deals with spin-wave modes in finite-size sam- 
ples of various geometries giving rise to ferromagnetic 
resonance (FMR). These modes were usually treated in 
the magnetostatic approximation, in which only the long- 
range dipole energy is retained in the dynamical equa- 
tions, while the exchange energy, as well as well as the 
retardation effects, are neglected 0,0- 

A special class of restricted-geometry spin modes is a 
surface magnetic wave, which propagates along the sur- 
face of a magnetically ordered material. In addition to a 
great physical interest of surface modes, they can be im- 
portant for various technological applications since 
for these waves the surface serves as a waveguide, which 
allows effective spin transport. The latter attracts now a 
lot of attention because of the prospects of spintronics. 

The most famous example of a magnetic surface wave 
is the Damon-Eshbach wave 2], which exists in slabs 
with the equilibrium magnetization parallel to the surface 
(parallel geometry). While the existence of this mode 
can be revealed within the magnetostatic approximation 
alone, it is not sufficient to obtain the correct spectrum. 
Thus, within this approximation, the surface mode lies 
above the bulk modes. It has been demonstrated Q that 
inclusion of the exchange energy can modify the results 
dramatically. Moreover, in the perpendicular geometry, 
a surface mode can exist only if the exchange energy is 
included. However, the analysis in Ref. Q| was only car- 
ried out numerically for some chosen parameter values. 
Ref. considered the effect of displacement currents (for 
the parallel geometry), with the conclusion that they also 
can be important, but the exchange energy was disre- 
garded in that work. 

This state of affairs, when it is clear that the magne- 
tostatic approximation is insufficient for the description 
of the surface wave spectrum, while no inclusion of the 



neglected parts has been done systematically, is part of 
the motivation for this work. Another goal of this work is 
to study the surface-wave spectrum for a superconduct- 
ing ferromagnet (SCFM) and to compare its behavior to 
that of an insulating ferromagnet and a nonmagnetic su- 
perconductor. SCFM's have attracted a lot of attention 
during the last decade in the context of unconventional 
superconductivity. The existence of spin modes can serve 
as a clear proof of magnetic order in these materials, 
where detection of ferromagnetism is hindered by screen- 
ing Meissner currents Q . Knowledge of the surface- wave 
properties, while interesting by itself, might be also use- 
ful for a design of new experimental techniques to study 
coexistence of ferromagnetism and superconductivity. 

This Letter studies the spectrum of surface magnetiza- 
tion waves on a semi-infinite medium in the perpendic- 
ular geometry, taking into account the dipole and spin- 
exchange energies and the retardation effects. Two types 
of materials are considered: SCFM and an insulating fer- 
romagnet. In the SCFM, (longitudinal) plasmons are ex- 
cited together with the transverse magnetic modes in the 
surface wave; however, their presence is only important 
for determination of the upper termination point of the 
branch. On the other hand, in the insulator, plasmons 
are absent, and the surface wave survives to much larger 
wave vectors. The retardation effects are important both 
for very small and very large wave vectors, the latter case 
being relevant only for the insulator, since in the SCFM 
the branch does not survive to this regime. The spectrum 
is calculated for both small and large wave vectors. For 
large wave vectors, the spectrum is the same for both 
types of ferromagnets, apart from the difference in the 
branch termination points. At small wave vectors, the 
spectra are different for the two types of materials. The 
surface branch for the SCFM survives down to zero wave 
vector and is separated from the bulk-mode continuum 
by a finite gap. The branch for the insulator is termi- 
nated at a very small wave vector, where it collides with 
the bulk-continuum bottom. 

We consider a semi-infinite uniformly magnetized 
medium with the spontaneous magnetization perpendic- 
ular to the surface. Physically this corresponds to a slab 
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FIG. 1: Propagation geometry 



whose thickness is large enough that the other surface 
can be neglected. Our goal will be to investigate surface 
waves in two cases: insulating medium vs. a supercon- 
ducting one. We assume no dissipation in both cases, 
since only under this assumption, as well as in the semi- 
infinite geometry, are surface waves well-defined. The 
setup of the system is shown in Fig. ^ The z axis to- 
gether with the spontaneous magnetization M are per- 
pendicular to the surface, while the x and y axes paral- 
lel to it. We also assume a magnetic anisotropy in the 
system of the easy-axis type which is strong enough in 
order to stabilize the system against the flip of the mag- 
netization into the x-y plane. We will be looking for 
surface waves propagating along the x direction with a 
wave vector ko- That is, we will be interested in solu- 
tions of the dynamical equations (specified below), made 
of plane- wave combinations ~ exp(— iut + ikir) such that 
k,j = fc x + qiZ, and ko is real, while qi are complex. 

The magnetization is governed by the Landau-Lifshitz 
equation, which gives for a frequency u> 0: 

-gM x m(a + j 2 k 2 ) + gM x b, (1) 

where m is the dynamical magnetization such that M 
\ I () + m (m is in the x-y plane); g is the gyromag- 
netic ratio; a the magnetic anisotropy constant (assum- 
ing a > 47r); 7 the exchange stiffness constant, and, 
finally, b is the magnetic induction excited by the os- 
cillating magnetization m. This induction is related to 
the magnetization by the generalized London equation, 
which is given for plane-wave solutions by 
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Here A is the superconducting penetration depth (for an 
insulator, A — > oo); K = lo/c the electromagnetic wave 
vector. The last term ~ K 2 in the denominator takes 
into account the displacement currents. The vector mi 
is the magnetization component transverse to the wave 
vector k: 

k 2 m± — fc 2 m — k(k- m) = q m x + k 2 m y — qkom x z. (3) 

Note that when the denominator in the RHS of Eq. @ 
vanishes, this equation does not specify anymore b, but, 



rather, fixes the direction of m, so that k 2 m± = 0. This 
situation happens for the Damon- Eshbach wave Q . How- 
ever, in our geometry this possibility can not be realized, 
since it would lead to a non-zero m z component. Sub- 
stituting the last two equations into Eq. (JIJ, we obtain a 
closed equation of motion for the magnetization: 
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where Q = Lo/gM and A~ 2 = A~ 2 — K 2 . The last term 
in the RHS breaks the rotational symmetry of the prob- 
lem and mixes different circular polarizations. From this 
equation, the spectrum of different modes composing the 
surface wave is given by 



n 2 = [5a + 7 2 fc 2 + 
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(4) 



This is a complicated self-consistent equation in which 
the frequency enters both the LHS and the RHS through 
the definition of A. For given f2 and ko, it gives four 
magnetic modes. 

In addition to these magnetic modes, there is also 
a plasma mode in the system, consisting of a longi- 
tudinal electric field alone. Its spectrum is given by 
lu 2 = oj 2 + c 2 k 2 where lo p = c/A and c p << c is the 
plasmon velocity originating from the electron gas com- 
pressibility. We will see that in the SCFM, plasmons are 
excited together with the magnetic modes in the surface 
wave. The electromagnetic (EM) field inside the sample 
is given by a combination of the above-mentioned five 
modes (or four for the insulator). The exact combina- 
tion is determined by coupling to the EM field outside 
the sample using the appropriate boundary conditions. 

The EM field outside the slab can be found using the 
Maxwell equations for the vacuum with the frequency 
lu and wave vector K = /c x + fc z z, where k z is nega- 
tive imaginary and, as before, K = lo/c. The obtained 
fields should be coupled to the fields inside the slab by 
requiring continuity of the EM fields parallel to the slab 
surface. For the insulator, continuity of the normal fields 
is then satisfied automatically, while for the SCFM, it 
is an independent condition which determines the ampli- 
tude of the plasmon mode. Eliminating the latter from 
the equations (for the SCFM), we obtain from these con- 
ditions the following relations for the magnetic modes: 
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where q p = y — A 2 c 2 /c 2 — fcg is the z component of the 
plasmon wave vector (for the insulator, the term with q p 
should be dropped as the plasmons are irrelevant). In 
addition to these, two more conditions are needed for 
the amplitudes of spin wave modes. A possible choice for 
these can be the condition of vanishing spin currents at 
the slab surface 0, which gives 



where An = fi — Q^. Its form is given in Fig. |5Ja). 

As fco — > 0, only the long- wavelength mode con- 
tributes to the surface wave, and the surface-mode spec- 
trum becomes il — Ck (l + Co/2), where Co is the sur- 
face impedance of SCFM in the zero-frequency limit [||: 
Co = —iKXy/Jl and 
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These four conditions together with the requirement that 
qi are positive imaginary and k z negative imaginary spec- 
ify a full system of equations whose solution yields the 
spectrum of the surface wave. 

We have solved these equations in the limit of small 
and large wave vector: fco << 7 _1 and fco >> 7" 1 and 
found important differences between the SCFM and the 
insulator in these two regimes. For a small wave vector, 
fco << 7 _1 , the plasmon contribution in the SCFM can 
be neglected as it is small by the parameter (c p /c) -1 / 2 . 
Moreover, in this limit, the spin-wave modes have a defi- 
nite circular polarization, and their dispersion to leading 
order is given by the usual dispersion of spin waves in 



This is just the surface impedance of a usual supercon- 
ductor with magnetic permeability fi, so the magnetic 
order plays no role in this regime except for adding some 
effective permeability. The corresponding mode, exist- 
ing on the surface of a conductor, is known as the Zen- 
neck wave |9j. If the medium is not a superconductor, 
but a metal, then the wave is damped. Note that the 
bottom of the conduction band, given at fco — > by 
fl m = 5a — 7 2 A~ 2 + 2-\/47r 7 A~ 1 , is separated from the 
surface-wave branch by a finite gap. 

For the insulator at fco << 7 -1 , the surface- wave fre- 
quency is Q ~ 5a. Around this frequency the modes have 
the following wave vectors: k\ 2 ~ 7 _1 fco; k^ ~ 7 , and 
&4 ~ K ~ l/C. To obtain the surface wave, it is enough 
to keep the first two modes and use them in the two equa- 



SCFM with k || z lal3: ft = ±[5a+j 2 q 2 +4Tr/(l+q 2 \ 2 )}. tions for my from Eqs , jr-p. The sur f ace -wave spectrum 



As long as the condition fco >> K holds (which is 
true except for very small wave vectors ~ l/C, where 
C = c/gMo is the light velocity in magnetic units), re- 
tardation effects can be neglected. Then Eqs. 15161 lead 
to 



is found from the condition qi = —q2, which is also the 
condition for the bottom of bulk mode continuum. Thus 
to leading order in 7fco, these two coincide and are given 
by 
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5a + Y / 87T7 2 (fcg-2[<5a/C*] 2 ). (11) 



In the leading order only the positively polarized modes 
(1,2) are excited, and the above equation becomes 



Qi+ 92+9i 92 + A~ 



- A 2 fc z 9i 92(91 +9 2 )/2, (7) 



from which a linear surface- wave dispersion is found: n 
flh + s I fco I, where 



(a- I )\/l + 7- 1 A % /Z^nT , 8) 

and fi h = ,5a-7 2 A- 2 /2 + 7A- 1 (4Tr + 7 2 A- 2 /4) 1 / 2 . This 
frequency is the frequency of (nonuniform) FMR for a 
thick SCFM slab ||. It is remarkable that the surface- 
wave spectrum contains information about it. 

At very small wave vectors, fco ~ K ~ l/C, retarda- 
tion is important. When it is taken into account, the 
spectrum gets modified, and instead of the linear ex- 
pression, it is given by fl — n^ + s(p z — K 2 /p z ) where 
p 2 z = k 2 — K 2 . Solving this, we obtain the SCFM surface- 
wave spectrum for fco << 7 -1 : 
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The limit C — > 00 and fco — > (while fc >> l/C) pro- 
duces n = 5a, which corresponds to FMR frequency for 
the insulator. 

Taking into account the next order, it is found that the 
surface wave is situated slightly below the bottom of bulk 
mode continuum, with a small separation n sw — nbuik = 
7 2 (i-L 2 — fcg) ps — 7 2 (fcp — 5a 2 /C 2 ). For very small wave 
vectors, fco ~ K, the retardation effects become impor- 
tant. At this region the branch is terminated, not surviv- 
ing to fco = 0. The reason for this is that the bulk- wave 
continuum at very small fco < K [Fig. Hfb)] is differ- 
ent from that for the superconductor [Fig. Ufa)], since 
photons with velocity c/^/Jm can propagate in an insu- 
lating ferromagnet, where fjn = a /5a. Namely, when 
fco = K 2 (l + a /5a) the mode 4 starts propagating. At 
this point, given by fco = fc^m = (a + Sa^^Sa 1 / 2 /C, 
the continuum bottom has a cusp, going down steeply for 
smaller fco, where it is no more given by q\ — —92, but by 
^4 = 0. The surface-wave branch collides with the con- 
tinuum bottom at k^ in and is terminated there. Note 
that near this point the penetration depth of the wave 
diverges, so the other surface might become important. 
The behavior of the branch for the insulator at small fco 
is shown in Fig. |2Jb). 
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FIG. 2: Surface wave dispersion f2(fco) (solid lines) at small 
fco (a) for SCFM; (b) for insulator. The shaded regions cor- 
respond to the bulk-mode continuum. The dashed line in (a) 
represents the linear approximation of the surface-wave spec- 
trum, when the retardation effects are neglected. 



Next we consider the behavior of the surface wave 
branch for large feo >> 7 ■ To leading order in 7^0, 
the spectrum can be found using the spin-wave approx- 
imation, in which only the exchange stiffness energy is 
retained, so that the frequency is given by fl — 7 2 fc 2 
(the spin- wave mode obtained in this approximation will 
be called mode 1). Then the surface- wave spectrum for 
both the insulator and SCFM is found from the condition 
qi = 0, leading to 



this point, given by k = k$g ax = 7- 
uum bottom collides with the surface-wave branch and 
terminates it. We see that for fc < the condition 

K < A -1 << k is satisfied, as was stated above. In the 
insulator, on the other hand, the surface wave branch sur- 
vives to much larger wave vectors. It is terminated when 
the mode 1 starts propagating, that is, when q\ = 0, 



_1 riy 2 , the contin- 
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= V2K = 2" 1 /2 7 - 2 C. 



(12) 



In conclusion, we have considered a surface wave in su- 
perconducting and insulating semi-infinite ferromagnets, 
taking into account the dipole and the exchange-stiffness 
energies, as well as the displacement currents. The spin- 
exchange stiffness is crucial for all regimes. The displace- 
ment currents are important at very small wave vectors 
for both types of materials and very large wave vectors, 
near the branch termination point, for the insulator. The 
presence of plasmons in the SCFM does not modify the 
spectrum, but cuts off the branch at the plasma fre- 
quency. For the insulator, the surf ace- wave spectrum 
at very small kg approaches the continuum band bot- 
tom, whereas in the SCFM, the branch remains sepa- 
rated from the continuum bottom by a finite gap (which 
depends on 7, A and Mo). Hence the wave in the lat- 
ter is more robust against smearing by dissipation. The 
magnetic order in the SCFM makes a profound influ- 
ence on the surface wave (except at very small wave vec- 
tors ko << K). Hence the surface wave spectrum might 
provide information about magnetic order in supercon- 
ductors, in particular, in unconventional superconductors 
with broken time-reversal symmetry. 

We acknowledge discussions with M. Golosovsky. This 
work has been supported by the grant of the Israel 
Academy of Sciences and Humanities. 



which is identical to the continuum bottom. In order 
to find a small separation between the continuum bot- 
tom and the surface wave, the full spin-wave dispersion 
relation, Eq. Q must be used. Then the wave vec- 
tor ki has a small imaginary component perpendicular 



to the surface, qi = iir-f 2 fc , so that Q s 
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These expressions are correct as long 



as fco >> K, which is satisfied for the SCFM, as dis- 
cussed below. For the insulator, this condition loses its 
validity near the (upper) termination point, and there 
qi = iirj- 2 (k 2 -2K 2 )/(k 2 -K 2 ) 3 / 2 so that n sw ~n bulk = 
-7r 2 7 - 2 (fc 2 - 2K 2 ) 2 /{k 2 - K 2 f. While the surface-wave 
dispersion is the same for the two considered types of ma- 
terials at large fcrj, the branches are terminated at points 
which are very different. Namely, for the SCFM, the 
branch termination point is determined by the condition 
that plasmons start propagating. This happens at the 
plasma frequency, when A~ 2 — > 0, that is, K 2 — > A~ 2 . At 



[1] A. I. Akhiezer, V. G. Bar'yakhtar, and S. V. Peletminskii, 
Spin Waves (North-Holland, Amsterdam, 1968). 

[2] R. W. Damon and J. R. Eshbach, J. Phys. Chem. Solids 
19, 308 (1960). 

[3] R. Bhandari and Y. Miyazaki, Jpn. J. Appl. Phys., 
Part 1 39, 3076 (2000); R. Bhandari and Y. Miyazaki, 
Jpn. J. Appl. Phys., Part 1 40, 3768 (2001). 

[4] R. E. De Wames and T. Wolfram, Appl. Phys. Lett. 
15, 297 (1969); R. E. De Wames and T. Wolfram, 
J. Appl. Phys. 41, 987 (1970). 

[5] A. Hartstein, E. Burstein, A. A. Maradudin, R. Brewer, 
and R. F. Wallis, J. Phys. C 6, 1266 (1973). 

[6] V. Braude and E. B. Sonin, Phys. Rev. Lett. 93, 117001 
(2004). 

[7] T. K. Ng and C. M. Varma, Phys. Rev. B 58, 11624 (1998). 
[8] V. Braude and E. B. Sonin, in preparation. 
[9] A. Sommerfeld, Partial Differential Equations in Physics 
(Academic Press, N.Y., 1964). 



